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Abstract 

Cosmic strings can be created in the early universe during symmetry-breaking phase 
transitions, such as might arise if the gauge structure of the standard model is extended 
by additional U{1) factors at high energies. Cosmic strings present in the early universe 
form a network of long horizon-length segments, as well as a population of closed string 
loops. The closed loops are unstable against decay, and can be a source of non-thermal 
particle production. In this work we compute the density of WIMP dark matter formed 
by the decay of gauge theory cosmic string loops derived from a network of long strings 
in the scaling regime or under the influence of frictional forces. We find that for 
symmetry breaking scales larger than lO^'^ GeV, this mechanism has the potential to 
account for the observed relic density of dark matter. For symmetry breaking scales 
lower than this, the density of dark matter created by loop decays from a scaling string 
network lies below the observed value. In particular, the cosmic strings originating 
from a U{1) gauge symmetry broken near the electroweak scale, that could lead to a 
massive Z' gauge boson observable at the LHC, produce a negligibly small dark matter 
relic density by this mechanism. 



1 Introduction 



Cosmic strings are one- dimensional topological defects [H [2|, [3] . They can be formed during 
symmetry breaking phase transitions in the early universe [1], such as might arise if the 
gauge symmetry group of the standard model (SM) is enlarged by additional U{1) factors 
at high energies [5], El [7]. Such gauge factors appear naturally in many models of physics 
beyond the SM. A typical cosmic string network consists of horizon-length long strings as 
well as smaller closed string loops. The long strings carry a net conserved topological charge 
and are stable, while closed string loops do not have a net charge and can decay away. In 
the present work we study the production of dark matter by the decays of cosmic string 
loops [a [9]. 

The evolution of cosmic string networks in the early universe has been studied extensively. 
It is found that the long strings in the network stretch with the Hubble expansion, as well 
as form closed string loops by the process of reconnection when they intersect themselves or 
each other. After an initial transient, these two processes balance out such that the string 
network reaches a scaling regime, in which the total energy density of the network makes 
up a small, constant fraction of the dominant matter or radiation energy density of the 
universe [TOl [TT| [T2| [T3l [TH [15] . This string energy fraction is largely independent of the 
initial conditions, and is about Gn relative to the critical density, where G = l/(87rMpj) is 
Newton's constant and fi is the string tension. In terms of the scale of spontaneous symmetry 
breaking rj from which the strings arise, the tension is on the order of 



Cosmic string scaling can be modified if the strings have significant interactions with the 
thermal background that lead to an effective frictional force on the strings. The importance 
of friction to the evolution of a string network depends on the pattern of symmetry breaking 
as well as the temperature of the surrounding plasma. 

A possible decay product of the string loops formed by a scaling cosmic string network 
is dark matter (DM) [16]. Indirect evidence for cold dark matter has been obtained from a 
number of sources. Together, they indicate a dark matter density of [17] 



This dark matter density can be explained by the existence of a new weakly-interacting 
stable particle with a mass on the order of the electroweak scale (WIMP) Such particles 
are common in extensions of the SM that stabilize the electroweak scale against quantum 
corrections such as supersymmetry [18], universal extra dimensions [19], and little Higgs 
models with T-parity |20j . 

String loops can produce DM and other particles in a number of ways. For local 
cosmic strings, corresponding to the spontaneous breakdown of a gauge symmetry, the direct 
emission of particles by the strings is thought to be suppressed [21] Instead, loops lose 
most of their energy by oscillating and emitting gravitational radiation [23] . However, when 

^However, see Refs. [131 [22] for arguments to the contrary. 



/i ~ rj'^. 



(1.1) 



nmih^ = 0.1143 ±0.0034 



(1.2) 
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the loop radius shrinks down to the order of the string width, it will self-annihilate into 
its constituent fields [8J. The subsequent decays of these states can produce dark matter. 
An even larger source of string annihilation and particle production by closed loops is the 
formation of cusps [24j. These are points on a loop where the string segment folds back on 
itself and briefly approaches the speed of light [25], [26]. In the vicinity of a cusp, a small 
portion of the string will self-annihilate creating particles [Ml [27] . In many simple solutions 
for the motion of a string loop, a cusp is generally found to occur about once per loop 
oscillation period [251 [26l [281 [211 [30l [31] . 

The amount of dark matter created by a cosmic string network depends on the initial size 
of the string loops that are formed by the network. While the evolution of the long horizon 
length strings is well understood, the details of loop formation are less clear. These details 
are closely related to the spectrum of small fluctuations on long strings. Significant advances 
have been made recently in this direction, both in numerical simulations [321 [331 [M] , as well 
as in analytic models Ref . [351 [33 [SZl [381 [SSI [IQ] . In the present work we will mostly adopt 
the results of Refs. [361 [23 [Ml [39] to characterize the initial loop size spectrum. 

In the picture of loop formation (in the scaling regime) that emerges from Refs. [361 [371 
[38l [39] , fluctuations on long strings are created near the horizon scale dn ^ t, which is also 
the scale that characterizes the long string network in the scaling regime. After they are 
formed, these fluctuations grow less quickly than the horizon, and thereby shrink relative to 
the characteristic scale of the long string network. The fluctuation spectrum that emerges 
is a power law in the fluctuation size that increases going to smaller scales. This power 
law is eventually cut off well below the horizon by gravitational radiation damping, which 
erases very small fluctuations. The cut-off occurs when the fluctuation size falls below the 
gravitational radiation scale [36] 

iGw = r{GfiY+'^t (1.3) 

where F ~ 50 is a constant, t is the cosmic time, and x = 0.10 during radiation and 0.25 
during the matter era. The small fluctuation spectrum on long strings is therefore peaked 
near low- This peak implies that lew sets the typical initial loop size ii [361 EZl [3H1 [39] . 

i^ ^ lew. (1.4) 

Both the number and the energy density of the loops formed are dominated by loops of this 
initial size|l However, recent simulations (and the analytic model of Ref. |40j) also point 
toward a significant loop population near the horizon scale, with ii ~ (0.1) tj [321 [331 [3^ - 
We will therefore consider larger loops as well in our analysis. 

The primary goal of this paper is to compute the dark matter density generated by a 
network of gauge theory cosmic strings in the scaling regime. We also calculate the dark 
matter generated when frictional interactions of the strings with the background plasma 
modify the evolution of the network. A necessary condition for the phenomenological 

^For smaller values of rj and at very early times, the gravitational damping length law can fall below the 
width of the string, w ^ r/^^. In this case, we expect that small fluctuations are cut off at £ ^ w > law by 
direct particle emission as suggested in Ref. [TH [22] ■ 
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viability of a cosmic string network is that the dark matter density it generates not exceed the 
observational bound. Our results therefore provide a constraint on theories of new physics 
that lead to the formation of cosmic strings, such as models containing new U{1) gauge 
groups that arise frequently in superstring theory constructions and in models of grand 
unification O [6l [7j. Our results also motivate certain classes of gauge extensions of the 
standard model that are able to generate the observed dark matter relic density by cosmic 
string loop decay. While we concentrate on the dark matter produced by a network of cosmic 
strings, our methods can also be applied to computing the densities of other cosmologically 
interesting particles arising from cosmic string loop decays, such as moduli and gravitinos. 

We find that for symmetry breaking scales rj larger than 10^° GeV, the decays of cosmic 
string loops derived from a network of long cosmic strings that are scaling or dominated 
by friction can potentially generate (more than) enough cold dark matter to account for 
the observed relic density. The precise density depends on the typical initial loop size and 
the effective branching fraction of the loop decays into cold dark matter, in addition to 
7]. Note that t] = 10^° GeV corresponds to G/i ~ 10~^^, for which the standard cosmic 
string signatures such as gravitational radiation and gravitational lensing are expected to be 
unobservably weak [21 [3]. When rj is much smaller than 10^° GeV, including those scales for 
which the massive Z' gauge boson associated with the symmetry breaking could be visible at 
the LHC, the density of non-thermal dark matter produced by this mechanism is well below 
the observed value. Thus, such models are not constrained by dark matter overproduction 
from scaling string loops. 

This paper is organized as follows. In Section [2] we derive a general formula for the 
production of dark matter by cosmic string loops. In Section [3], we apply this formula to 
compute the amount of dark matter produced by string loops in the scaling regime. In 
Section H] we extend our results to compute the dark matter density from cosmic string 
loops when friction plays an important role in the evolution of the string network. Finally, 
Section [5] is reserved for our conclusions. 

Before proceeding, let us point out that the production of dark matter by the decays 
of cosmic string loops was considered previously for ordinary Abelian Higgs model cosmic 
strings in Ref. [8], for cosmic strings derived from a supersymmetric flat direction in Ref. |41j . 
and for cosmic strings associated with aspects of superstring theory in Refs. [9l HJl US] . We 
concentrate on ordinary cosmic strings in the present work, and we update and extend the 
results of Ref. [8] in a number of ways. Most importantly, we focus on cusping as the primary 
source of particle production by cosmic strings. We also make use of the recent results on 
the size distribution of string loops when they are formed from Refs. |36l [371 [Ml [39] , and 
we consider the evolution of cosmic string networks both with and without friction. When 
friction is relevant, we extend Ref. [8] by using the analytic model of Ref. [13] to describe 
the long string network. 
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2 A Formula for Dark Matter from Cosmic Strings 



Let us denote the initial invariant length of a string loop formed at cosmic time tj by ii. 
The invariant length £ of a loop is defined in relation to its energy in the cosmological frame 
by 

Eioop = l^i. (2.1) 

If the loop is boosted with speed u, the invariant length £ will exceed the proper length of 
the loop in its rest frame by a factor of 7 = 1/ ^1 — u'^. 

The key quantity characterizing loop formation is 

r{ei,ti)de,du, (2.2) 

the number density of string loops formed in the time interval (tj, ti + dti) with initial length 
in the range + dii). The form of the function r{ii,ti) is constrained by the evolution 
of the long string network. In particular, the total rate at which energy is transferred from 
the long string network to loops is typically a known quantity [13]. We can use this fact to 
impose the constraint 



dpi 



oop 



diipiir{ii,t,), (2.3) 



dti 

where dpioop/ dt is the rate of energy transfer into loops from the long string network. 

Consider the evolution of a loop of initial size ii formed at time ti. The length of this 
loop is described by the function 

iit;ii,ti) (2.4) 

which is the solution to the equation 

di 

Pj^=-Ptou with i{ti-ii,ti) =1,, (2.5) 

where Ptot is the total rate of energy loss from the loop. We will assume that Ptot is positive 
so that £ is monotonically decreasing in time. Under this assumption, the loop will eventually 
decay away completely at the time tco{ii,ti), defined implicitly by 

= £{t,o;ii,ti). (2.6) 

While the loop is decaying, it will emit a fraction of its energy in the form of (cold) dark 
matter. We will write this fraction as Pdm < Ptot- 

To find the total rate of dark matter production by the collection of string loops, consider 
loops with initial size in the range {ii, ii + dii) formed in the time interval (tj, ti + dti). There 
are r{ii,ti) diidti such loops per unit volume initially. As time goes on, this collection of 
loops is diluted by the cosmic expansion by a factor of a~^. This is the only modification of 
their density until they decay away completely at time tco{ii,ti). At time i, this collection 



5 



of loops will produce dark matter at the rate Pdm- The dark matter will thermalize if it is 
produced before the freeze-out time i < tjo, and redshift as after this time. Thus, we 
have 

- /, NT 3 r /rx n 3 

a[ti 



Ap 



DM 



Mi 



to 



dti 



tco i^^iiti) 



dt r{ii,ti 



fo 



DM 



a(t) 
a{to) 



Qitco-tfo). (2.7) 



The integral over ti runs between the string network formation time and the present time 
to ^ 4 X 10^1 GeV-\ Typically, t„ is on the order of 



f ~ 



M, 



PI 



(2i 



where rj is the scale of spontaneous symmetry breaking. The integration over the decay time 
t runs from the freeze-out time tfo to either t^o or t^, whichever is smaller. Thus, we have 
defined tco{ii,ti) by 



to 



tco < to 
tco > to 



(2.9) 



Any restriction on the integration limits of ii are encoded in the support of the function 
r{ii,ti). 

It is often the case that Pdm and Ptot depend on i only through i = i(i;ii,ti). If so, it 
is convenient to change the integration variable from t to i. The corresponding Jacobian is 
simply {d£/di)~^ = —p/Ptot, where Ptot is the total power released by a loop of length i. It 
follows that 



DM 



to 



dti 



diriii,ti 

(^i ) 



ajti) 
a{to) 



1 3 



P, 



DM 



P 



tot 



The limits on the i integration depend on the functions 



as well as 



^ foi^^i^ ti 



(^j ; ti 



^(tfoj^iyti^ . ti < tjo 



ti > t 



fo 



: t^o{ii,ti) <to 

^{to'Ai^ti) : tco{ii,ti) > to 



(2.10) 



(2.11) 



(2.12) 



Eqs. (12. 7p and (12.101) are our main results. We will apply them to two interesting special 
cases below. 



3 Dark Matter Production from a Scaling Network 

As a first application of our main result, Eq. (12. 7p . we compute the dark matter density pro- 
duced by a network of cosmic strings in the scaling regime due to the cusping of closed string 
loops. To simplify the analysis, we focus on monochromatic loop formation distributions with 

ii = ati. (3.1) 
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This relation implies r{£i,ti) oc S{£i — aU). More general distributions can be obtained by 
introducing a weight function and summing the final result over different values of a. 

In the string scaling regime, the rate at which the long string network transfers energy 
into loops is [13] 

= C/itr', (3.2) 



-3 

dU ' 

with C ~ 10 a constant characterizing the mean properties of the long string network. 
Imposing the constraint of Eq. (I2.3p on the loop formation rate r{ii,ti), we obtain 

r{i„U) = ^t;U{ii-ati). (3.3) 

a 

All that remains to do is to specify the evolution of the loop length and the power emitted 
as dark matter, and apply Eq. (12.71) . 

Cosmic string loops lose energy to gravitational radiation as well as cusping, and shrink 
as a result. The rate of energy emission into gravity waves is j3] 

Pgrav = TGfX^ (3.4) 

where /i is the string tension, G = 1/ (STrMpJ, and F ~ 50 is a dimensionless constant 



Loops also lose energy when they form cusps, which are points on the loops that briefiy 
fold back upon themselves and approach the speed of hght. Near the apex of a cusp, a 
portion of the string overlaps itself and self-annihilates. Summing over many cusps, the net 
rate at which a loop loses energy to cusping is given by [27] 

rW 

Pcusp = fJ'Pcd J, (3.5) 

where w ~ f]^^ is the string width, and Pc is the probability for a cusp to form per period 
of loop oscillation. Several studies suggest Pc — 1 t26l [28l [29l [30l [3T|. Comparing the cusp 
power of Eq. (13.51) to the gravitational wave power of Eq. (13. 4p . we see that cusping is the 
dominant energy-loss mechanism by loops when they are shorter than £ < with 

(3.6) 

The particles produced by cusping can decay into dark matter. Provided Pc ^ TGfi, we 
expect cusping to be the dominant source of particle production and dark matter from 
cosmic strings. 

If loops lose energy to gravitational radiation and cusping, the loop length evolves 
according to 

= -^tot = -TGfi^ - f^PcJj- (3.7) 
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The solution of this equation subject to the initial condition i{ti) = ii is given implicitly by 



t-U 



- 2 




+ 2 In 



(3.8) 



Suppose that a fraction ecusp of the energy emitted by a string loop cusp (eventually) 
takes the form of cold dark matter. It follows that 



DM ^cusp Pcusp- 



(3.9) 



Putting this into Eq. fl2.10p and making use of Eq. (13.31) . we find the dark matter density 
due to cusping to be 



DM 



^cusp 



Pcuspi^^) 



Mi / dU / dir{£i,U) { — ] fi 

Jtr, Jije^,h) \«o/ Ptot{f^) 



(3.10) 



to 



dti — t, 



a 



ao 



2 (^cusp 




In (l + 



e.fo(oiti,ti) 



The B functions in this expression account for the fact that only those loops that decay 
after tfo and before to can contribute to the dark matter density. Numerically, we find that 
the integrand of Eq. fl3.10p is a steeply falling function of ti that is cut off at small ti by 
the Q{C.fo) condition. It follows that the dominant contribution to the dark matter density 
comes from loops formed at the smallest possible value of ti such that they decay shortly 
after tfo. In this sense, the dark matter density is created nearly instantaneously at tfo- 



To be concrete, we evaluate Eq. (13.101) assuming a weakly-interacting massive dark matter 
particle with a freeze-out time of t/o = 2 x 10^^ GeV~^. This corresponds approximately 
to a freeze-out temperature of 5 GeV, which is in the range expected for a stable, weakly- 
interacting particle with a mass on the order of 100 GeV. We take the universe to be radiation 
dominated prior to tfo-, and we assume a standard cosmological evolution afterward. We 
also set the loop network parameter to C = 10, the cusping probability to = 1, and the 
branching fraction into dark matter equal to unity, ecusp = 1- In realistic models ecusp can 
be much smaller than unity. Our results can therefore be interpreted as providing an upper 
bound on ecusp within the underlying gauge theory model. 

In Fig. [1] we show the dark matter relic density due to cusping as a function of the initial 
loop size parameter a, normalized to FG/i, with the other parameter values as described 
above. We normalize a to FG/x because, in the absence of cusping, a loop is long-lived relative 
to the Hubble time at formation (~ ti) provided a/TGfi > 1, and short-lived otherwise. From 
this plot, we see that increasing the initial loop size a well above FG/i increases the final dark 
matter density. The reason for this is that the integration over ti in Eq. (I3.10p is dominated 
by the earhest times for which loops produced at ti decay after tfo- Since larger loops are 
longer-lived, loops contributing to the dark matter density can be formed at earlier times 
when the integrand of Eq. (I3.10p is larger. As a/TGfi decreases, the dark matter density 
curves flatten out for lower values of the VEV rj. In these flat regions, cusping dominates the 
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Figure 1: Dark matter density due to loop cusping for ecusp = 1, Pc = 1, C = 10, and 
tfo = 2 X 10^^ GeV~^ as a function of the initial loop size parameter a. The various lines 
correspond to different values of the symmetry breaking VEV rj. 



evolution of the loops contributing to the dark matter, and these loops are short-lived even 
though a/TGfi > 1. When the VEV is larger, t] > 10^^ GeV, gravitational radiation becomes 
more important than cusping. Loops then become short-lived only when a/T Gfi < 1. When 
a becomes very small, cusping again takes over and the curves flatten out once more. 

In Fig. [2] we show the dependence of the dark matter density from string cusping as 
a function of the symmetry breaking VEV t]. The other string model parameters are as 
described above. This plot shows a general increase in the dark matter density up to large 
values of t], and then a fall-off. Short-lived loops, with a < F G/i, and smaller values of 
7] decay mostly through cusping at time tfo. Using Eq. (12.101) . it is possible to show that 
Qdm increases as 77^ in this case. As r] is increased further, gravitational radiation becomes 
more important than cusping, and loops lose most of their energy to gravity waves instead 
of dark matter. Thus, the dark matter density falls for these very large values of 77, with 
^DM oc (ar^)"^/^ in this regime. The precise value of r] at which the cross-over occurs 
depends on the value of a, with smaller initial loops being more prone to cusping. For very 
large initial loops, a = 0.1, a similar cross-over occurs. 

In general, we expect the loop cusping probability pc to be on the order of unityjfl 
However, it is also interesting to look at smaller values of pc, such as might arise on large 
string loops with many kinks [SQl [31]. Fig. [3] shows the effect of reducing pc for small 
(a = FGyu) and large {a = 0.1) initial loop sizes, and a range of values of the symmetry 

•^The backreaction from string annihilation at a cusp can suppress its reoccurrence, but this need not 
prevent the formation of new cusps |27j . 
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Figure 2: Dark matter density due to loop cusping for e^usp = l5C = 10, Pc = l, and t/o = 
2x 10^^ GeV^^ as a function of the symmetry breaking VEV r/. The various curves correspond 
to different values of the initial loop size parameter a. The black dashed horizontal line 
denotes the observed dark matter relic density. 



breaking VEV r^. All other parameters are as in the previous plots. For a = F G/i the loops 
are necessarily short-lived. In this case, cusping is the dominant energy loss mechanism at 
time tjo for i] < 10^^ GeV, and remains so even for lower values of Pc- Thus, reducing 
does not alter the resulting dark matter density. When ri is larger, gravitational radiation 
dominates the energy loss at time t/o, and reducing -pc therefore decreases the fraction of 
loop energy released as dark matter. With a = 0.1 string loops can be long-lived. For 
smaller values of < 10^^ GeV, lowering can enhance the dark matter density. This 
occurs because cusping is the dominant energy loss mechanism around time tj^ for these 
loops. Reducing thus further increases the lifetime of these long-lived loops which, as 
discussed above, enhances the final dark matter density. For larger values of 77, gravitational 
radiation dominates at time t/o? so that reducing only decreases the fraction of energy 
released as dark matter without significantly altering the loop lifetime. 

We have only computed here the contribution to the dark matter density from loops 
formed in the scaling regime. Loops created in the string-forming phase transition at time t,, 
and while the string network was evolving towards scaling will give an additional contribution 
to the dark matter density if they decay after tjo- The approach of a string network to scaling 
depends on the details of the phase transition, and a study of this process is beyond the scope 
of the present work. However, we can derive the condition for loops formed in the phase 
transition to decay before t/o. We obtain 
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Figure 3: Dark matter density due to loop cusping for ecusp = 1, C = 10, and tfo = 2 x 
10^^ GeV^^ as a function of the cusp formation probability pc for several values of the 
symmetry breaking scale t]. The panel on the left corresponds to small initial loops with 
a = FG/i, "while the panel on the right corresponds to large initial loops with a = 0.1. The 
black dashed horizontal line indicates the observed dark matter relic density. 

If this bound is satisfied and scaling is attained quickly after the network is formed at t^, 
loops formed in the scaling regime will give the dominant contribution to the dark matter 
density. This bound also indicates that our results are applicable to the phenomenologically 
interesting case of a U{1) gauge symmetry broken only slightly above the electroweak scale 
provided ~ 1 and a -C 0.1. 

Before moving on, let us mention two additional effects that we have not yet taken into 
account. The first is that we have treated the loops as having no net motion relative to 
the thermal background, even though it is found in Refs. [361 EH ES [39j that small loops 
(a < F G/i) can be highly boosted. This boost does not alter the cusping power, although 
the expression for this power in Eq. (13. 5p picks up a factor of 7^/^ (for short-lived loops), 
where 7 = — u"^ is the boost factor. The decay products will also be boosted. Unless 

the loops are very tiny, however, this boost is not expected to greatly modify our results. 

A second and potentially much larger effect that we have not yet included is that the 
decay products from a loop cusp can be very highly boosted, even if the loop is initially at 
rest. The size of this boost 7 was estimated in Ref. [27] to be 




(3.12) 



On account of this boost, the cusp decay products including the dark matter might not 
thermalize immediately upon production. If not, they may thermalize at a later time, or 
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even remain as hot dark matter. 

To estimate the effect of this possibly very large boost, we adopt a simple picture of the 
cusp decay products as consisting of a bunch of light states (including WIMP dark matter) 
moving together with a boost 7 relative to the thermal background. In the rest frame of the 
bunch we assume very little relative motion, and that a fraction ei of the total bunch energy 
consists of dark matter particles of mass m ~ 100 GeV. The dark matter in the bunch will 
thermalize when aritot ^ H, where ntot ~ is the total plasma number density of light 
states and a is the cross-section for the dark matter to scatter off the thermal background. 
To estimate the thermalization time, we take this cross-section to be 

_ 1 1 _ 1 1 

8tt s 8tt '-y'^im? 

With this cross-section, the dark matter decay products from a cusp on a small loop with 
a = TGfi formed at time tjo will reach kinetic equilibrium nearly instantaneously after the 
cusp occurs when the symmetry breaking scale lies below 77 < 10^^ GeV. In this case, our 
previous estimates for the dark matter density remain valid. 

If the cusp decay products do not thermalize instantaneously, they may still thermalize at 
a later time tx as their momentum redshifts. To compute the final density of cold dark matter 
in this case, it is only necessary to add a factor of a(t)/a(tr) to the integrand of Eq. (12.71) 
and modify the limits of integration such that tx > tfo rather than t > tfo- Doing so, we find 
that the resulting dark matter density coincides with what we found previously (without 
the boost) up to an overall factor of order unity even when instantaneous thermalization 
does not occur. The additional redshift a(t)/a(tr) is cancelled nearly exactly by the delay 
in termalization, which allows loops contributing to the dark matter density to be formed at 
earlier times when the rate of loop production was larger. Similar to before, we also find that 
the dominant contribution to the cold dark matter density comes from loop decay products 
that thermalize shortly after tfo- 

The prospect of highly boosted decay products also raises several new issues. Since a large 
fraction of the energy of the cusp typically goes into the kinetic energy of the decay products, 
the fraction of this energy that ends up as cold dark matter can be suppressed. If kinetic 
thermalization proceeds exclusively through elastic collisions, we would have ecusp = ei/7(^r), 
where 7(tT) is the boost factor at the time of thermalization. The actual value of ecusp can be 
larger than this if inelastic collisions during thermalization generate additional dark matter 
particles, but in general we expect a significant suppression of ecusp relative to the non- 
boosted case. Highly boosted loop decay products can also generate relics that have not 
thermalized by the present time (or the time of matter-radiation equality). In general, the 
resulting densities of hot dark matter are suppressed by factors of Gfi and are safely small 
for lower values of the symmetry breaking scale. However, a reliable precise estimate is 
complicated by the fact that for extremely large boosts, the cusp decay products will naively 
have super-Planckian momenta. A full analysis of this issue is beyond the scope of the 
present work. 

In summary, our results indicate that the density of WIMP cold dark matter generated by 
string loops from cusping with pc = 1 in the scaling regime is safely below the observed value 
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(3.13) 



for 1] < 10^° GeV {G^ ~ 10^^^). For values of the symmetry-breaking scale larger than this, 
the string-induced dark matter density depends on the initial loop size. Loops of initial size 
a = T (G/i)^-^, motivated by the results of Refs. [361 EH [SHJ [39], can create as much or more 
than the observed dark matter density for f] ^ 5 x 10^^ GeV (G/x ~ 10~^^) if the branching 
fraction ecusp is on the order of unity. Larger loops, with a ~ 0.1, can lead to a much greater 
dark matter density. If such large loops are typical, the branching fraction into dark matter 
Ecusp must be significantly less than unity if the cosmic strings are to avoid overproducing 
dark matter. This is natural for large loops, whose cusps generate very highly boosted decay 
products, leading to much of the cusp energy being lost to kinetic energy that is transferred 
to the thermal bath rather than cold dark matter. Let us also emphasize that for models 
with a new U{1) gauge symmetry broken only slightly above the electroweak scale, the cold 
dark matter produced by the corresponding cosmic string network in the scaling regime is 
negligibly small for both large and small initial loop sizes. 

4 Dark Matter Production with Friction 

The evolution of a cosmic string network can be modified if the strings have significant 
interactions with the thermal background. As we will discuss below, the relevance of these 
interactions to the string network depends on the details of the symmetry breaking from 
which the cosmic strings arise. Such interactions, when present, tend to slow the motion of 
the strings by creating an effective frictional force on them H5l [13] . This in turn changes 
the density and rate of growth of the network. The frictional forces on strings decrease as the 
universe cools, and eventually become unimportant relevant to the Hubble damping from 
the expansion of spacetime. Frictional effects also change the way cosmic string loops form 
and decay, and can enhance the total rate of loop formation. In this section we apply the 
result of Eq. (12.71) to compute the density of dark matter created by a cosmic string network 
evolving under the influence of frictional forces. 

For local (gauge) cosmic strings, the dominant interaction between the strings and the 
thermal background comes from Aharonov-Bohm scattering [IB]. This results from the 
phase change experienced by the charged particle as it is transported around the string. The 
effective frictional force induced by this scattering can be characterized by a friction length 
£f. Friction becomes unimportant when if{t) grows larger than the Hubble length. The 
friction length due to Aharonov-Bohm scattering is given by [HI, [T5] 



where the dimensionless quantity (3 is [Bj 



with the sum running over relativistic degrees of freedom, and ba = 1(3/4) for bosons 
(fermions). The value of Ua is related to the charge Qa of the light particle species a. If 




(4.2) 



a 
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the underlying U{1) gauge symmetry (subgroup) is broken by the condensation of a field 
with charge Q, it is equal to Qa/Q- Therefore we expect f3 to be of order unity when some 
of the i>a are non-integer, and zero otherwise. When Aharonov-Bohm scattering vanishes, 
Everett scattering of charged particles off the strings will be the dominant source of frictional 
interactions [17|. The corresponding friction length is similar to that for (non-trivial) 
Aharonov-Bohm scattering but is enhanced by a factor of In^^T/i]). Frictional interactions 
will be completely irrelevant when all the light states in the theory are uncharged under the 
broken gauge group. In the present section, we will assume that Aharonov-Bohm scattering 
is the dominant source of frictional interactions with (3 = 1. 

Frictional effects on the long string network decouple when the friction length grows 
larger than the Hubble length. This occurs at the time defined by the relation 

H{t,) = 1/ifit,). (4.3) 

Friction is only relevant for the long string network when t < t^. For (3 = 0{1) and radiation 
domination, t^ has the parametric siz^ 

(. ^ ^, (4.4) 

where t] denotes the symmetry breaking VEV. This is parametrically larger than the typical 
(radiation-era) formation time t^ ~ M^i/rj^. 

The evolution of a cosmic string network in the presence of friction was studied in 
Ref. In their analytic model, the long string network is characterized by an effective 

correlation length L and a mean velocity v. The energy density of the network is given in 
terms of these variables as 

Poo = (4.5) 

If the initial string density is larger than the scaling density, as would be expected if 
the symmetry breaking phase transition is second-order or weakly first-order, the string 
network evolves very quickly to the Kibble regime [H |T3]. In this regime, with the universe 
assumed to be radiation dominated, the string network variables L and v have the parametric 
dependences 

Lit) ^ (^-j t (4.6) 

Ht) - (-) . (4.7) 

The Kibble regime only lasts while t < t^. From Eqs. (14.1 1) and (14.41) we see that the friction 
length grows as 

] t. (4.8) 



^ Due to the many uncertainties involved in the description of cosmic strings within the friction-dominated 
regime, we only list and use here the leading parametric dependences of the string network parameters. 
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As t approaches t*, the friction length catches up to the long string length L as well as the 
horizon, and the string network transitions into the usual scaling regime with L oc t and 
1/ ~ 1. 

During the friction-dominated Kibble regime, the rate at which energy is transferred from 
the long string network to loops is on the order of 

This is parametrically larger than during the scaling regime, as can be seen by comparing 
with Eq. fl3.2p . Once a loop is formed, its subsequent evolution in the Kibble regime is 
also considerably different than during scaling. In particular, the typical initial loop size as 
well as the subsequent loop evolution are both strongly modified. Both of these effects can 
modify the resulting dark matter density. 



4.1 Loop Production and Evolution with Friction 

The evolution of linear perturbations on long strings and closed string loops in the presence 
of friction was studied in Refs. |^H5] . In Ref. [15] it was found that linear fluctuations on a 
long string of wavelength larger than if are overdamped and stretched. For wavelengths much 
smaller than if, the damping time is on the order oiif, which is much longer than the typical 
period of oscillation but much shorter than the Hubble time. Hence these small fluctuations 
oscillate and lose energy to friction very quickly relative to the Hubble time. Given these 
results and the picture of loop formation of Ref. [39j, we expect that the typical initial loop 
size during friction is on the order of Fluctuations smaller than if are damped out 
quickly, while those larger than if grow more slowly than the long string correlation length 
L, and therefore shrink relative to L. Thus, we expect that fluctuations build up near if, 
which in turn sets the typical size of a loop when it is formed. 



To model the evolution of string loops during friction, we will take the results of Ref. 
for the evolution of a circular loop to be representative of the evolution of general loops. 
(Indeed, friction tends to make the loops more circular.) Ref. pS] finds that loops smaller 
than if oscillate essentially freely and lose their energy to the thermal background according 
to 



^-/x-. (4.10) 

friction f 

It follows that such loops lose energy over the time scale if, much less than the Hubble 
time in the friction regime. We do not expect these interactions between the strings and the 
thermal background to be a significant source of dark matter. 

The motion of loops larger than ^/ is overdamped. They evolve according to 

RR^ — ^ 4.11 
2a 



^We emphasize however that the picture of loop formation obtained in Ref. [551 '^^^ developed under the 
assumption of long string scaling, and did not consider the effects of friction. 
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where R is the comoving coordinate radius of the loop (and aR corresponds to the physical 
loop radius). The solution of Eq. (14.111) implies that loops of initial size smaller than the 
long string scale L shrink down to size if in less than about a Hubble time. Once they do, 
the overdamping approximation of Eq. (14. lip breaks down, and the loops begin to oscillate 
and decay away. 

In summary, loops are formed in the friction regime with a typical initial length close to 
if. The loop formation rate per unit volume per unit length is 

r(£„t,)~ (^1^ t-U{i,~if{U)), U<U. (4.12) 

Any loop formed with initial length larger than if will shrink down to length in less than 
about a Hubble time. Loops smaller than if execute underdamped oscillations, transferring 
their energy to the thermal background and decaying away over the time scale if. The loop 
length i in this regime evolves according to 

- = ---TG,^p.^-. (4.13) 

Here, the first term comes from friction, the second from gravitational radiation, and the 
third from cusping. 

Before going on to compute the dark matter density generated by the decaying loops, 
let us make note of the fact that at the end of the friction-dominated Kibble regime, the 
long string network is smoothed out nearly all the way up to the Hubble scale. In the loop 
formation picture of Ref. (33, ESI ES], small fluctuations on long strings giving rise to 
loops originate from fluctuations of Hubble size that have slowly shrunk. Thus it will take 
some time for small scale structure to build up on the long strings, and the typical initial loop 
size will be initially larger than the gravitational damping length law- We follow Ref. pj] 
and model this transitional period by writing the initial typical loop size parameter as 

2 + ait/U)^ 

with the exponent ^ ~ 1. A nai've application of the results of Refs. EH EEl EH] suggests 
that ^ = 0.9 in the radiation era. We will study a range of values of ^. 



4.2 Dark Matter with Friction 

The discussion above provides all the ingredients needed to evaluate the dark matter density 
created by a string network in the presence of friction using Eq. (12.71) . We do so here under 
the assumption that loop cusping is the dominant source of dark matter from the strings. Our 
results are presented in Fig. IH where we show the dark matter density due to cosmic strings 
as a function of the symmetry breaking VEV rj. In making this plot, we have also set C = 10, 
Pc = 1, Gcusp = 1, and we have again taken the freeze-out time to be t/o = 2 x 10^^ GeV~^. 
We have also fixed an overall prefactor in Eq. (14.121) (equal to 2^/(2-f-a)) to ensure that the 
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Figure 4: Dark matter density due to loop cusping for ecusp = 1, Pc = 1, and t/o = 2 x 
10^^ GeV~^ as a function of the symmetry breaking VEV r] in the presence of friction. The 
various curves correspond to different values of the initial loop size parameter a. The black 
dashed horizontal line indicates the observed dark matter relic density. 



loop formation rate is continuous as ti crosses t*. The different curves in Fig. H] correspond 
to different values of the scaling regime loop size parameter a and the exponent ^ appearing 
in Eq. dHU). 

The dark matter density curves in Fig. |4] show three distinct regions, with the transitions 
between them occurring around t] = 10^ GeV and r] = 10^^ GeV. For values of t] well above 
10^^ GeV, the curves coincide with those obtained in the absence of friction. Such large 
values of t] imply a value of t^, that is much smaller than the freeze-out time tfo, so that all 
loops generated while friction is relevant decay away before tfo- On the other hand, when t] 
lies below 10^^ GeV, the most important contribution to the dark matter comes from loops 
formed while friction dominates the network evolution, ti < t*. 

In the region 10^ GeV ^ f] ^ 10^^ GeV most of the dark matter is produced by loops 
that are formed in the friction era, with tj < t*. As a result, the dark matter density is 
independent of the scaling value of a for this range of rj. The largest contribution to the 
dark matter in this region comes from loops that are also long-lived, with tj < tfo. This 
enhances the amount of dark matter formed because, with the loop distribution function of 
Eq. (14.121) . the integrand of the ti integral in Eq. (12.71) is a rapidly decreasing function of 
ti. Increasing t] in this region leads to larger initial loop sizes that further extend the loop 
lifetime. However, the minimal value of ti for which loops decay after tfo decreases more 
slowly with increasing rj than t^,. The sharp transition near t] = 10^^ GeV occurs when these 
two quantities become equal. When this happens, increasing t] further decreases the initial 
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loop size (see Eq. fl4.14p ) for the dominant DM loops, and these loops go quickly from being 
long-lived to being short-lived, with the dominant loops being formed near ti = tfo > t^. 
The transition is more gradual when the exponent ^ in Eq. (14.141) is less than ^ = 1.0, which 
can be seen by comparing with the curve for ^ = 0.7. 

For 7] < 10^ GeV, the initial loop size if becomes small enough that the loops are short- 
lived, decaying away within a Hubble time. Thus, the majority of the dark matter produced 
for these smaller values of rj comes from loops formed near the freeze-out time, ti ~ tfo, with 
initial size close to if (tfo). As rj decreases below 10^ GeV, we find that a larger fraction of 
the energy of each loop is lost to cusping, thereby increasing the dark matter density. Even 
so, the dark matter relic density generated by cosmic strings for 77 ~ 10^ GeV is a negligibly 
small fraction of the observed value. 



5 Conclusions 

We have investigated the cold dark matter density created by the decays of loops of gauge 
cosmic strings. Dark matter is produced by string loops when they form cusps. At a 
cusp, a small portion of the string loop annihilates into its constituent fields, which can 
then cascade down to lighter states such as dark matter particles. Our results provide 
constraints on extensions of the gauge symmetry group of the standard model that give rise 
to cosmic strings in the early universe. The string loops that give decay to dark matter 
are themselves created continually by the network of long, horizon-length strings. We have 
studied the amount of dark matter generated when this network is in the scaling regime, 
and when its evolution is dominated by frictional forces. Both cases are physically relevant 
as the presence or absence of significant frictional interactions depends on the details of the 
symmetry breaking from which the strings originate. 

The dark matter density generated by a string network in the scaling regime can po- 
tentially be enough to explain the observed relic density when the breaking VEV rj exceeds 
10^° GeV. The amount of dark matter produced by this mechanism is greatest when the 
initial loop size approaches a significant fraction of the cosmological horizon, although much 
smaller initial loop sizes can also generate sizeable densities of dark matter. However, larger 
loops have cusp decay products that are highly boosted, and this can reduce the fraction of 
the loop energy that becomes cold dark matter. For values of rj below 10^° GeV, the amount 
cold dark matter created by a scaling string network is always well below the observed value. 

String networks that are strongly influenced by frictional interactions can also generate 
a dark matter density that is equal to or greater than the observed value. This can occur 
if the symmetry breaking scale is greater than about 10^^ GeV. For very large values of 
the symmetry breaking scale, above about 10^^ GeV, frictional effects decouple well before 
the non-thermal dark matter is created; the frictional interactions do not alter the resulting 
dark matter density. When the symmetry breaking scale is much below 10^^ GeV, the 
contributions to the dark matter density in the presence of friction are very small. This 
contrasts with and is less constraining than the result of Ref. jSj, in which the effects of 
cusping on the decays of loops were not included in computing the dark matter density. 
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We conclude from our investigation that the non-thermal cold dark matter relic density is 
safely below the observed value when the symmetry breaking scale is less than rj ^ 10^'' GeV. 
For larger symmetry breaking scales, the non-thermal dark matter density produced by 
decaying string loops can potentially be as large or larger than the observed value. Thus, 
our results also put constraints on new gauge symmetries broken at very high scales based 
on the requirement that they do not overproduce dark matter. These constraints are 
model-dependent, in that they depend on the effective branching fraction of the decaying 
string loops into dark matter, but they can be more severe than the constraints from the 
more traditional cosmic string signatures such as gravitational radiation and gravitational 
lensing [21 [3]. For gauge symmetries broken at lower scales that could be probed by the LHC, 
we find that the contribution to the dark matter density from their corresponding scaling 
string loops is negligibly small. 
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